A step toward CNO solar neutrinos detection in liquid scintillators 
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Abstract 



^ The detection of CNO solar neutrinos in ultrapure liquid scintillator detectors is limited by the back- 

ground produced by Bismuth-210 nuclei that undergo /3-decay to Polonium-210 with a lifetime of ~ 7 days. 
Polonium-210 nuclei are unstable and decay with a lifetime equal to ~ 200 days emitting a particles that 
can be also detected. In this letter, we show that the Bi-210 background can be determined by looking at 

the time evolution of a— decay rate of Po-210, provided that a particle detection efficiency is stable over the 

(-H data acquisition period and external sources of Po-210 are negligible. A sufficient accuracy can be obtained 

^ I in a relatively short time. As an example, if the initial Po-210 event rate is ~ 2000 cpd/lOOton or lower, a 

Ji. Borexino-like detector could start discerning CNO neutrino signal from Bi-210 background in At ~ lyr. 

> 
in 
m 
cn 

^' 

o 

T — I 

X 



1 



1 Introduction 



One of the main goals of the present and next generation ultrapure Hquid scintillator detectors, such 
as KamLAND [1], Borexino [2], SNO+ [3J and LENA is the determination of the neutrino fluxes 
produced by the CNO cycle in the Sun. 

Despite being sub-dominant in the Sun, the CNO cycle p] has a key role in astrophysics, being 
the prominent source of energy in more massive stars and in advanced evolutionary stages of solar like 
stars, see [6|. The evaluation of CNO efficiency is connected with various interesting problems, like 
e.g. the determination of globular clusters age [7] from which we extract a lower limit to the age of the 
Universe. At the moment, we still miss a direct observational evidence for CNO energy generation in 
the Sun. We only have a loose upper limit on CNO luminosity, LcnOi obtained by combining the results 
of the various solar neutrino experiments. The most recent Solar Standard Model (SSM) [8] predicts 
the CNO contribution to the Sun's luminosity to be equal to about 0.7%. Experimental bounds report 
-^^CNO < 3% at 99% C.L., see e.g. |9l [lOl [11]. The detection of CNO solar neutrinos would clearly 
provide a direct test of the CNO cycle efficiency. 

The measurement of the CNO solar neutrino flux can also provide clues to solve the so called "solar 
composition problem". The flux is, in fact, directly related to the abundance of carbon, nitrogen and 
oxygen in the Sun. The photospheric abundances of these (and other) heavy elements have been recently 
re-determined \T2\ [T3l I14j . indicating that the sun metallicity is lower than previously assumed ^15j . 
Solar models that incorporate these lower abundances are no more able to reproduce the helioseismic 
results. Detailed studies have been done to solve this discrepancy, see e.g. [iHl [IT] , but a definitive 
solution of this problem is still missing. 

The above points show the importance of CNO solar neutrinos detection which, however, is a very 
difficult task. Not only the flux is relatively low, but also their energy is not large. The neutrinos 
produced in the CNO cycle have continuous energy spectra with endpoints at about ~ 1.5 MeV[^ 
Differently from the monochromatic Be and pep solar neutrinos, they do not produce specific spectral 
features that permit to extract them unambiguously from the background event spectrum in high 
purity liquid scintillators. In particular, as it is shown in Fig{T| the electrons produced by the /3-decay 
of Bismuth-210 to Polonium-210 have a spectrum that is similar to that produced by CNO neutrinos. 
As a consequence, spectral fits are able to determine only combined "Bismuth+CNO" contribution, as 
it is done e.g. by Borexino in [21 111] . 

In order to remove this degeneracy, we propose a simple method to determine the Bi-210 decay 
rate which is based on the relationship between the Bi-210 and Po-210 abundances. Polonium-210, 
which is the Bismuth-210 daughter, is unstable and decays with a lifetime rpo ~ 200 days emitting a 
monochromatic a particle that can be easily detected. In the absence of Bismuth-210, the a— decay 
rate of Po-210 nuclei follows the exponential decay law, npo(t) oc exp(— t/rpo). The deviations from this 
behaviour can be used to determine the /3— decay rate of Bi-210 nuclei. The only required assumptions 
are that the a particle detection efficiency is stable and external sources of Po-210 are negligible during 
the data acquisition period. 

The plan of the paper is the following. In the next section, we determine the relationship between 
the a— decay rate of Po-210, npo(t), and the /3— decay rate of Bi-210, nBi(t). In sectjsj we estimate the 
accuracy AnBi of the Bi-210 decay rate determination as a function of the detector mass M, the data 

^ The dominant components of the CNO neutrino fiux are the so called N-neutrinos produced by ^^N — >^'' C + e+ + !/c 
with an endpoint at E,y = 1.20 MeV and the 0-neutrinos produced by ^^O — N + e+ + !/c with an endpoint at Ei, — 1.73 
MeV. 
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Figure 1: The expected event spectrum in solar neutrino liquid scintillator detectors calculated by assuming: 
the solar neutrino fluxes predicted by fBli the oscillation parameters corresponding to the LMA-MSW flavour 
oscillation solution liyj/ : background levels and detector energy resolution comparable to those obtained by Borexino 
detector Jllf . See sect^for details. 



acquisition period At and the initial polonium decay rate npo,o- In sectjlj we discuss the implications 
of the proposed approach for CNO neutrino signal extraction. In sectjsj we summarize our results. 

2 The relationship between Polonium and Bismuth 

Bismuth-210 and Polonium-210 are both daughter of ^^^U. Bismuth-210 is produced by the slow decay 
of Lead-210 which has a lifetime equal to rpb = 32.3 y. It then undergoes a f3~ decay to Polonium-210: 

with a lifetime tbi = 7.232 d. The electrons produced in the decay have a continuous spectrum with an 
endpoint -Emax = 1.16 MeV that lays over the event spectrum produced by CNO neutrinos, see Fig{T] 
In the analysis of the 192 days data |llj . Borexino obtained an event rate equal to about 20 cpd/lOOtons 
[llj for the combined contribution provided by Bismuth+CNO neutrinos. For reference, we consider 
that the expected CNO neutrino signal is equal to about 5 cpd/lOOton as predicted by the SSM and 
the LMA-MSW scenario (see sect J4] for details). 

Polonium-210 is also unstable and undergoes a decay to Lead-206: 

210po^206p^^^ (2) 

with a lifetime rpo = 199.634 d. The a particles produced in the decay have an energy Ea = 5.3 MeV, 
but they are observed at a much lower effective energy, due to the large quenching factor of alpha 
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particles in liquid scintillators. They produce a sharp peak superimposed to the low energy part of the 
''Be neutrino event rate at a visible energy Svis — 0.5 MeV, see Fig{l} 

Naively, one could expect that secular equilibrium between Bismuth-210 and Polonium-210 is ob- 
tained, which would imply that the rate of /3 decays of Bi-210 is equal to the rate of a decays of Po-210. 
However, this is not the case in practical situations. Borexino, in fact, observes a Po-210 abundance 
which is about a factor 100 larger than what implied by secular equilibrium and that corresponds to 
an initial alpha decay rate of the order ~ 8 x 10^ cpd/lOOton [18j. The large counting rate permits to 
determine the Polonium activity with a high accuracy. 

We remark that the a particles emitted in the decay, beside producing the very distinctive spectral 
feature shown in Fig|l| can be directly identified by using a pulse shape discrimination. As it is shown 
in I2i lUj a high purity liquid scintillator could have very good a — /3 discrimination properties. A 
figure-of-merit based on the Gatti discrimination technique [20] has been used in Borexino to produce 
an a-subtracted spectrum where the Po-210 background is removed. 

The relationship between the Po-210 and the Bi-210 abundances can be quantified in a very simple 
way. We have that: 

— t;— = \ V ^Po(rj Ki) 

at Tpo TBi 

where N-p^ify (-^Bi(i)) is the number of Polonium-210 (Bismuth-210) nuclei per unit mass in the detector 
at a time t, and S'po(t) indicates any possible external source of polonium-210 (i.e. not related to ^^^U 
decay chain). The solution of Eq. ([3]) can be written as: 

A^Po(t) = iVpo,o exp(-t/rpo) + rpo /^^^ + 5po(t)\ (4) 

\ TBI I 

where A^po,o is the number of Po-210 nuclei per unit mass at the time t = and the symbol {f{t)) 
indicates the time "average": 

(fit)) = — f dt' fit - t') exp(-t7rpo) (5) 

TPo JO 

for the generic function f{t). 

We assume that, during a data acquisition period, after liquid scintillator purification, external 
sources of Polonium can be neglected, i.e. S-poif) ~ 0. In this assumption, the a— decay rate of 
Polonium and the /3— decay rate of Bismuth, given by: 

npo(t) = Npo{t)/Tpo 

riBiit) = NBi{t)/TBi (6) 

respectively, follow the relation: 

npo(i) = npofi exp(-t/Tpo) + (nBi(t)) . (7) 

The above equation provides all the ingredients that are necessary for our analysis. It shows that, if 
we are able to measure rapo(t) as a function of time with a high accuracy, we can provide an estimate 
for nBi(i), which is the main background for CNO neutrinos detection. 
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Figure 2: The la expected accuracy Anei in Bi-210 event rate determination. 

3 Determination of the Bi-210 event rate 



The main observation of this letter is that, due to the high statistics of Po-210 events, the beta-activity 
of Bi-210 nuclei can be extracted from the time evolution of Polonium-210 counting rate in a relatively 
short time (i.e. comparable with Polonium lifetime). It is not necessary, in particular, to wait that 
secular equilibrium between Polonium-210 and Bismuth-210 is reached, which would require much 
longer times. 

This point can be quantified with straightforward analytical considerations by making the following 
exercise. We assume that the Bi-210 event rate is constanl|^ obtaining from Eq.Q: 

n-po{t) = [npo,o - nBi] exp(-t/Tpo) + nei . (8) 

We imagine to collect the Po-210 event^ for a total period equal to At. We then divide this time 
interval into two large bins from [0, At/2] and [At/2, At]. According to the time evolution implied by 
Eq.®, the number of events collected in the two bins are: 



iVi 

N2 



eM 
eM 



At \ 

(npo,o - riBi) rpo ( 1 - e ^-po + nei — 



{npofi - nBi) Tpo e ^^po ( 1 - e ^^po I + ubi 



At 



(9) 
(10) 



where M indicates the detector mass and e indicates the detection efficiency (averaged over the inte- 
gration period) that we assume not to vary significantly. 



This hypothesis corresponds to assuming that there are no sources of Lead-210 and Bismuth-210 in the detector and 
it is done here only to consider a specific situation and to illustrate quantitatively the proposed method. We remark that, 
in the analysis of the experimental data, it is not necessary to make a-prion assumptions on the time evolution of nBi{t) 
since this can be determined from the experimental data themselves, as it is discussed in sect j4] and appendix. 

•^The Po-210 counting rate can be determined either by fitting the low energy part of the event spectrum or by selecting 
Po-210 events by pulse shape discrimination tecnique, see [20]. The optimal approach has to be choosen on experimental 
basis. 
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Table 1: The absolute uncertainty Atibi in the Bi-210 event rate determination, reconstructed by performing a 
fit over an observation period At to simulated data. All rates are expressed in cpd/lOOton. We assume that the 
true Bi-210 event rate is nei = 20cpd/100ton. 



In the absence of the Bismuth-210 contribution, the ratio of events collected in the two bins is equal 
to N2/N1 = exp (— At/2Tpo). The deviations of N2/N1 from this value can be used to measure the beta 
activity of Bi-210. We have, in fact: 



N2- Nie 



e MnBi 



At 



(11) 



By propagating the statistical errors = V^2 and AA'^i = \/Ni, one is able to estimate the accuracy 
AriBi of the determination of the Bi-210 decay rate. In the assumption that npoo ^ '^Bii we obtain: 



AnBi 



"-Po,0 



/(At) 



rpoM 

where we considered e ~ 1. The function /(At) is explicitly given as: 



(12) 



/(At) 



2tpo 
At 



l + e 



At 
"2To 



1 



At 



(13) 



We see that the relevant parameters are the initial Po-210 event rate npo,0) the detector mass M and the 
total time of the measure At. In Figj2] we show the iso-countour lines for AriBi in the plane (npo,0) At) 
for two values of the detector mass M = 100 ton and M = 1000 ton. For an initial Polonium activity 
equal to npo,o = 2000 cpd/100 ton, we reach an accuracy comparable to the expected CNO signal, 
ncNO - 5 cpd/100 ton, after a time At ~ 300d (At ~ 150d) for a 100 ton (1000 ton) detector. We 
remark that an important assumption in the above analysis is that the detection efficiency e is known 
and stable over the period of the measure. Variations and/or errors in e clearly provide additional 
sources of uncertainties. 

These analytic estimates have been checked numerically by simulating random Po-210 events, 
according to the time evolution expressed by Eq.Q, in detectors with masses M = 100 ton and 
M = 1000 ton. We assumed that the Bi-210 event rate is equal to riBi = 20 cpd/lOOton, while 
the initial Po-210 event rate have been varied in the range npo,o = 2000 — 8000 cpd/lOOton in order to 
study the dependence of our final results on the initial polonium contamination. We binned the data 
over five days and then extracted ubi and npo,o by fitting the data over an observation period At with 
the functional form in Eq.Q. The accuracy AriBi for the reconstructed Bi-210 event rate is reported 
m tab{lj We see that the numerical results agree within ~ 15% with the analytical expression (|12l). 
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4 Implications for CNO solar neutrino extraction 



The possibility to determine the Bi-210 event rate represents an important step forward in the direction 
of CNO solar neutrino detection, as can be immediately understood by looking at the expected event 
spectrum in liquid scintillator detectors in FigjlJ The spectrum was calculated by assuming the solar 
neutrino fluxes predicted by [8j and the oscillation parameters corresponding to the LMA-MSW solution 
of the solar neutrino problem given by [H]. The integrated rates are n-Qe = 50cpd/100ton from 
beryllium neutrinos, Up^p = 2.7cpd/100ton for pep neutrinos and = 0.48 cpd/lOOton for boron 
neutrinos. The CNO neutrino signal is taken as ncNO = 5.1 cpd/lOOton that corresponds to the 
prediction obtained by using the high photospheric metal abundances of [15j. For comparison, by using 
the low photospheric values of [13] one obtains the lower rate ncNO = 3.6 cpd/lOOton. The background 
levels have been estimated considering the Borexino results as a reference (see e.g. [llj). The Kripton-85 
and Carbon-11 event rate are chosen as nxr = 25cpd/100ton and ncii = 25 cpd/lOOton, respectively. 
The Bi-210 event rate is taken as ubi = 20cpd/100ton, while the (initial) Po-210 event rate is chosen 
as npo,o = 2000cpd/100ton. This value is lower than what initially obtained in Borexino |18| but we 
believe that it represents a reasonable goal for next runs and next generation experiments. Finally, 
we modelled the detector energy resolution by a Gaussian function with an energy dependent width 
a{E) = 0.05 MeV MeV tUj. 

We see that the CNO neutrino fluxes are not expected to produce significant and/or recognizable 
features anywhere in spectrum (see Fig{T]). As a consequence, spectral fits are unable to constrain the 
CNO signal. This is emphasized in the left panel of Figj3| where we display the results of a spectral fit to 
simulated data in the plane (nsi, ?^CNo)- In the simulation, we considered a detector mass M = 100 ton 
and a total observation period equal to At = 1 yr, during which all the components have been assumed 
to be constant in time except the Po-210 contribution that was modulated in time according to Eq.Q. 
In the fit, we considered all the signals and backgrounds as free parameters, except the pp, pep and 
Boron neutrinos contributions which have been fixed to their expected values. We see from Fig jS] that 
the CNO signal is basically unconstrained. The origin of the strong anti-correlation with the Bi-210 
background is understood on a quantitative basis by considering that the Bi-210 and CNO components 
give the dominant contribution to the spectrum only in a narrow energy range at ii^ ~ 0.8 MeV. The 
signal produced in this window is proportional to the combination: 

= ncNO + ^, mi (14) 

where the factor ^ = 0.67 is the ratio between the normalized Bi-210 and CNO spectra evaluated at 
E ~ 0.8 MeV. The data allow, thus, a good determination of with a poor reconstruction of ncNO 
and riBi- 

The degeneracy between Bi-210 and CNO neutrinos is removed when we use the time evolution of 
the Po-210 contribution to determine the Bi-210 background. By performing a fit to the simulated data 
in the domain of energy and time and taking advantage of Eq.Q between Po-210 and Bi-210 rates, we 
obtain the considerable improvement in the CNO signal determination that is shown in the right panel 



of Fig c The final uncertainty is AncNO — 2.1 cpd/lOOton and it is essentially determined by the 
accuracy of the Bi-210 event rate determination according to AncNO ~ C AnBi, as can be estimated 
by using the Anei values given in tabjl] and/or in Eq.(12). 



*This plot has been obtained by fitting simultaneously the time and energy distribution of simulated data without 
considering the possibility to identify Po-210 events by pulse shape discrimination [20], since this additional information 
do not introduce relevant improvements in our analysis. This possibility, however, should not be overlooked because it 
may became important in the analysis of real data, depending on the experimental conditions. 
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The above results indicate that present generation hquid scintihator experiment, hke e.g. Borexino, 
already have the potential to probe the CNO neutrino flux, provided that they are stable for sufficiently 
long time (~ lyr) and/or the initial polonium contamination can be made sufficiently low. Considering, 



moreover, that the uncertainty scales as AnBi oc \J npo^o/M, we note that: i) future Kton-scale detectors, 
like e.g. SNO+, will be able to start discriminating between high and low metallicity solar models; 
ii) the initial Po-210 contamination npo^o ~ 8000 cpd/lOOton reported Borexino by p!8] corresponds, 
anyhow, to la extraction of CNO solar neutrinos. 

A final remark is important. In our analysis, we assumed that the Bi-210 event rate is constant in 
time. This assumption was done to consider a specific situation but it is not necessary in the analysis of 
real data, since the time evolution of nBi(i) can be determined from the experimental data themselves. 
A simple strategy can be the following. By performing spectral fits to the signal at different times t, 



we determine the quantity n/3(i) defined in Eq.(14). This allow us to obtain Bi-210 rate as a function 
of time, according to: 

miit) = - [npit) - ncNo] (15) 
? 

where the unknown factor ncNO is constant in time (we neglect the small seasonal modulation of 
neutrino fluxes due to earth orbit eccentricity). The presence of a time varying nBi(i) does not invalidate 
the possibility to perform our analysis, since the only requisite for the validity of Eq.Q is that there 
are no external Polonium-210 sources. We can use, in fact, Eq.(15) in Eq.([7]) to obtain the following 
expression: 

npo(t) - (n^(t)) = [npofi + ncNo] exp(-t/Tpo) - ncNO (16) 

that allows to determine ncNO from a two parameters fit to the "observable" quantity npo(t) — {np{t)) 
with no additional assumptions. Alternatively, we can compare the behaviour of np{t) with expecta- 
tions. A constant value is expected when external sources of Bismuth-210 and Lead-210 are negligible. 
A linear growth of nBi(i) (and thus of np{t)) is, instead, produced by a slow increase of the Radon 
contamination of the detector (see appendix for details). We can assume an arbitrary functional form 
for n-Q\{t) that agrees with the observed np{t) and fit the data in the time and energy domain, using 
Eq.Q in place of Eq.Q to discriminate Bi-210 from CNO neutrino contribution. 

5 Conclusions 

In this letter we have proposed a simple and general method that allows to measure the flux of CNO 
neutrinos in massive high radio-purity liquid scintillator detectors by means of determining the Bi-210 
background. We summarize here the main points of our analysis: 

i) We have discussed the relationship between Polonium-210 and Bismuth-210, see Eq.Q, showing 
that the deviations of the a— decay rate of Po-210 from the normal exponential decay law, npo(t) oc 
exp(— t/rpo), can be used to determine the beta activity of Bi-210. 

ii) We have estimated the expected accuracy of the Bi-210 event rate determination as a function of the 
relevant parameters that are the initial polonium contamination of the detector, the detector mass and 
the total observational time. Our results are summarized by Eq.(12) and/or tabjlj A good accuracy 



can be obtained in relatively short time, comparable to the Po-210 lifetime rpo ~ 200 d. 
Hi) We have discussed the implications of our approach for CNO solar neutrino detection. We have 
shown that, if the initial Po-210 event rate is ~ 2000 cpd/lOOton or lower, a Borexino-like detector 
could start discerning CNO neutrino signal in Ai ~ lyr. Future Kton-scale detectors, like e.g. SNO-I-, 
in the same time interval, could begin to discriminate between high and low metallicity solar models. 
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^Bi (cpd/lOOton) ^Bi (cpd/lOOton) 



Figure 3: The l-a and 2-a allowed regions obtained by a fit to simulated data, assuming a detector mass 
M = 100 ton and an observation period At = lyr. In the left panel, we consider only information contained in 
the energy distribution of the events. In the right panel, we perform a fit in the domain of energy and time, taking 
into account the relationship between Po-210 and Bi-210 expressed by Eq.^. The blue dashed line corresponds 
to the condition np = const, see Eq.(14-). 



The determination of the Bi-210 event rate could represent an important step toward CNO solar 
neutrinos detection in high purity liquid scintillator experiments. The only required assumptions are 
that a— particle detection efficiency is stable over the data acquisitions period and external sources of 
Po-210 are negligible. A special effort on the experimental side in this direction is worthwhile in the 
next future. 
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A Lead-Bismuth-Polonium 210 decay chain 

The decay chain is described by the system of differential equations: 

^ = -« + ^ + 52(t) (17) 

where, to keep notation simple, we use the index 1 is for Lead-210, the index 2 for Bismuth-210 and 
the index 3 for Polonium-210. The functions Si{t) describe source terms and the initial conditions are 
expressed as: 

Ni{0) = Ni,o (18) 
These equations can be treated as a decoupled system of equations in the form 

Njjt) ^ Njjt) 
dt Ti 

that is readily solved as follows: 



+m (19) 



Ni{t) = Nie n +n{m)i (20) 

where we used the notation 

m)i = - fe-'-^mdl^ (21) 

Ti Jo 

that can be thought of as a suitable average over an interval of time of the order of Ti preceding the 
time t. Some explicit expressions are given in the following. 

The solution of the homogeneous system (i.e., with source terms set to zero) with the boundary 
condition iVj(O) = A^^^o is: 

Ni{t)= e'^Ni^o 



£21 
-<:21 



e '^3 

where we used the notation 



(22) 



(23) 



and considered the inequality n » T3 » r2. The inhomogeneous (full) system admits the following 
solution satisfying the boundary condition (5iVj(0) = 0: 

5Ni{t) = Ti{Si{t))i 

6N2{t) = T2{S2{t) + {Si{t))i)2 (24) 
SNsit) = MSsit) + {S2it) + (Si(t))i)2)3 

Thus, due to the linearity of the system, the most general solution in this case is just 

Ni{t) = Ni{t) + 5Ni{t) (25) 
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As an application, we discuss the case when we do not have sources of Bismuth or Polonium but 
only of Lead (this corresponds to the practical situation in which the Radon content of the detector 
increases with time, since Radon-222 converts rapidly in Lead-210). We consider the simplest case, 
namely when the source is a constant: 



Slit) = Si, S2{t) = S^it) = 



(26) 



In this assumption, eqs.(24) can be explicitly calculated obtaining: 
6Ni{t) = Smil - e^n) 

6N2it) = SiT2 (l 



e ^1 h e ^2 



(27) 



6N3it) = Sirs 



1 



e ^1 



1 

(l-e2l)(l-€3l) 



e ^2 



£21623 



(l-<:2l)(l-e23) 



+ e ^3 



£31 



(l-£23)(l-£3l) 



Now consider the fact T2 <^ <^ ti and the limit T2 <^ t <^ ti, which applies to the typical conditions 
of data taking, t being on the scale of one year. We get the approximated solutions: 



Afpb(t) «iVpb,o + i(5pb 
iVBi(t) « ^ iVpb(t) 



-^Pb,0 ' 
Tpb ' 



(28) 



iVpo(t) « iVpcO exp (-^) + rpo 



\ •^Bi /po 



where we used the notations adopted in the paper to increase readability. If the Lead-210 source is 
sufficiently intense, i.e. Sp^ > -^'pb,o/''"Pb) we expect a linear increase of the Bi-210 beta activity in 
the detector, nBi(i) = ^Bi(i)/''"Bi- When the Lead-210 source can be neglected, we obtain instead 
nBi(i) — const. 
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